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METHOD AND APPARATUS FOR EXPONENTIALLY CONVERGENT 
THERAPY EFFECTIVENES MONITORING USING DNA 
MICROARRAY BASED VIRAL LOAD MEASUREMENTS 
Field of the Invention 

The invention generally relates to techniques for monitoring the effectiveness of 
medical therapies and dosage formulations, and in particular to techniques for monitoring 
therapy effectiveness using viral load measurements. 



10 Background of the Invention 

It is often desirable to determine the effectiveness of therapies, such as those directed 
against viral infections, including therapies involving individual drugs, combinations of 
drugs, or other related therapies. One conventional technique for monitoring the 

1 5 effectiveness of a viral infection therapy is to measure and track a viral load associated with 
the viral infection, wherein the viral load is a measurement of a number of copies of the virus 
within a given quantity of blood, such per milliliter of blood. The therapy is deemed 
effective if the viral load is decreased as a result of the therapy. A determination of whether 
any particular therapy is effective is helpful in determining the appropriate therapy for a 

20 particular patient and also for determining whether a particular therapy is effective for an 
entire class of patients. The latter is typically necessary in order to obtain FDA approval of 
any new drug or medical device therapy. Viral load monitoring is also useful for research 
purposes such as for assessing the effectiveness of new antiviral compounds determine, for 
example, whether it is useful to continue developing particular antiviral compounds or to 

25 attempt to gain FDA market approval. 



A test to determine the viral load can be done with blood drawn from T-cells or from 
other standard sources. The viral load is typically reported either as an absolute number, i.e., 
the number of virus particles per milliliter of blood or on a logarithmic scale. Likewise, 
decreases in viral load are reported in absolute numbers, logarithmic scales, or as 
percentages. 

It should be noted though that a viral load captures only a fraction of the total virus in 
the body of the patient, i.e., it tracks only the quantity of circulating virus. However, viral 
load is an important clinical marker because the quantity of circulating virus is the most 
important factor in determining disease outcome, as changes in the viral load occur prior to 
changes in other detectable factors, such as CD4 levels. Indeed, a measurement of the viral 
load is rapidly becoming the acceptable method for predicting clinical progression of certain 
diseases such as HIV. 

Insofar as HIV is concerned, HIV-progression studies have indicated a significant 
correlation between the risk of acquiring AIDS and an initial HIV baseline viral load level. 
In addition to predicting the risk of disease progression, viral load testing is useful in 
predicting the risk of transmission. In this regard, infected individuals with higher viral load 
are more likely to transmit the virus than others. 

Currently, there are several different systems for monitoring viral load including 
quantitive polymerase chain reaction (PCR) and nucleic acid hybridization. Herein, the term 
viral load refers to any virological measurement using RNA, DNA, or p24 antigen in plasma. 
Note that viral RNA is a more sensitive marker than p24 antigen. p24 antigen has been 
shown to be detectable in less than 50% asymptomatic individuals. Moreover, levels of 
viral RNA rise and fall more rapidly than levels of CD4+ lymphocytes. Hence, changes in 
infection can be detected more quickly using viral load studies based upon viral RNA than 
using CD4 studies. Moreover, viral load values have to date proven to be an earlier and 
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better predictor of long term patient outcome than CD4-cell counts. Thus, viral load 
determinations are rapidly becoming an important decision aid for anti-retro viral therapy 
and disease management. Viral load studies, however, have not yet completely replaced 
CD+ analysis in part because viral load only monitors the progress of the virus during 
5 infection whereas CD4+ analysis monitors the immune system directly. Nevertheless, even 
where CD4+ analysis is effective, viral load measurements can supplement information 
provided by the CD4 counts. For example, an individual undergoing long term treatment 
may appear stable based upon the observation of clinical parameters and CD4 counts. 
However, the viral load of the patient may nevertheless be increasing. Hence, a 
10 measurement of the viral load can potentially assist a physician in determining whether to 
change therapy despite the appearance of long term stability based upon CD4 counts. 

Thus, viral load measurements are very useful. However, there remains considerable 
room for improvement. One problem with current viral load measurements is that the 
threshold level for detection, i.e., the nadir of detection, is about 400-500 copies per 
1 5 milliliter. Hence, currently, if the viral load is below 400-500 copies per milliliter, the virus 
is undetectable. The virus may nevertheless remain within the body. Indeed, considerable 
quantities of the virus may remain within the lymph system. Accordingly, it would be 
desirable to provide an improved method for measuring viral load which permits viral load 
levels of less than 400-500 copies per milliliter to be reliably detected. 
20 Another problem with current viral load measurement techniques is that the 

techniques are typically only effective for detecting exponential changes in viral loads. In 
other words, current techniques will only reliably detect circumstances wherein the viral load 
increases or decreases by an order of magnitude, such by a factor of 10. In other cases, viral 
load measurements only detect a difference between undetectable levels of the virus and 
25 detectable levels of the virus. As can be appreciated, it would be highly desirable to provide 
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an improved method for tracking changes in viral load which does not require an exponential 
change in the viral load for detection or which does not require a change from an 
undetectable level to a detectable level. Indeed, with current techniques, an exponential or 
sub-exponential change in the viral load results only in a linear change in the parameters used 
5 to measure the viral load. It would instead be highly desirable to provide a method for 
monitoring the viral load which converts a linear change in the viral load into an exponential 
change within the parameters being measured to thereby permit very slight variations in viral 
load to be reliably detected. In other words, current viral load detection techniques are useful 
only as a qualitative estimator, rather than as a quantative estimator. 

10 One reason that current viral load measurements do not reliably track small scale 

fluctuations in the actual number of viruses is that a significant uncertainty in the 
measurements often occurs. As a result, individual viral load measurements have little 
statistical significance and a relatively large number of measurements must be made before 
any statistically significant conclusions can be drawn. As can be appreciated it would be 

1 5 desirable to provide a viral load detection technique which can reliably measure the viral 
load such that the statistical error associated with a single viral load measurement is 
relatively low to permit individual viral load measurements to be more effectively exploited. 

Moreover, because individual viral load measurements are not particularly significant 
when using current methods, treatment decisions for individual patients based upon the viral 

20 load measurements must be based only upon long term changes or trends in the viral load 
resulting in a delay in any decision to change therapy. It would be highly desirable to 
provide an improved method for measuring and tracking viral load such that treatment 
decisions can be made much more quickly based upon short term trends of measured viral 
load. 

25 As noted above, the current nadir of viral load detectability is at 400-500 copies of 



the virus per milliliter. Anything below that level is deemed to be undetectable. Currently 
the most successful and potent multi-drug therapies are able to suppress viral load below that 
level of detection in about 80-90 percent of patients. Thereafter, viral load is no longer an 
effective indicator of therapy. By providing a viral load monitoring technique which reduces 
the nadir of detectability significantly, the relative effectiveness of different multi-drug 
therapies can be more effectively compared. Indeed, new FDA guidelines for providing 
accelerated approval of a new drug containing regimen requires that the regimen suppress the 
viral load below the current nadir of detection in about 80 to 85 percent of cases. If the new 
regimen suppresses the viral load to undetectable levels in less than 80 to 85 percent of the 
cases, the new drug will gain accelerated approval only if it has other redeeming qualities 
such as a preferable dosing regimen (such as only once or twice per day), a favorable side 
effect profile, or a favorable resistance or cross-resistance profile. Thus, the ability of a 
regimen to suppress the viral load below the level of detection is an important factor in FDA 
approval. However, because the level of detectability remains relatively high, full approval 
is currently not granted by the FDA solely based upon the ability of the regimen to suppress 
the viral load below the minimum level of detection. Rather, for full approval, the FDA may 
require a further demonstration of the durability of the regimen, i.e., a demonstration that the 
drug regimen suppresses the viral load below the level of detectability and keeps it below the 
level of detectability for some period of time. 

As can be appreciated, if a new viral load measurement and tracking technique were 
developed which could reliably detect viral load at levels much lower than the current nadir 
of detectability, the FDA may be able, using the new technique, to much more precisely 
determine the effectiveness of a drug regimen for the purposes of granting approval such that 
a demonstration of the redeeming evalities will no longer be necessary. 

For all of these reasons, it would be highly desirable to provide an improved 



technique for measuring and tracking viral load capable of providing much more precise and 
reliable estimates of the viral load and in particular capable of reducing the nadir of 
detectability significantly. The present invention is directed to this end. 

Summary of the Invention 

In accordance with a first aspect of the invention, a method is provided for 
determining the effectiveness of a therapy, such as an anti-viral therapy, by analyzing biochip 
output patterns generated from biological samples taken at different sampling times from a 
patient undergoing the therapy. In accordance with the method, a viral diffusion curve 
associated with a therapy of interest is generated and each of the output patterns 
representative of hybridization activity is then mapped to coordinates on the viral diffusion 
curve using fractal filtering. A degree of convergence between the mapped coordinates on the 
viral diffusion curve is determined. Then, a determination is made as to whether the therapy 
of interest has been effective based upon the degree of convergence. 

In an exemplary embodiment, the viral diffusion curve is spatially parameterized such 
that samples map to coordinates near the curve maxima, if the viral load is increasing (i.e., 
therapy or dosage is ineffective). In this manner, any correlation between rate and extent of 
convergence across different patient samples is exploited to provide a quantitative and 
qualitative estimate of therapy effectiveness. 

Also in the exemplary embodiment, the biological sample is a DNA sample. The 
output pattern of the biochip is quantized as a dot spectrogram. The viral diffusion curve is 
generated by inputting parameters representative of viral load studies for the therapy of 
interest, generating a preliminary viral diffusion curve based upon the viral load studies; and 
then calibrating a degree of directional causality in the preliminary viral diffusion curve to 
yield the viral diffusion curve. The parameters representative of the viral load studies 
include one or more of baseline viral load (BVL) set point measurements at which detection 



is achieved, B VL at which therapy is recommended and viral load markers at which dosage 
therapy is recommended. The step of generating the preliminary viral diffusion curve is 
performed by selecting a canonical equation representative of the viral diffusion curve, 
determining expectation and mean response parameters for use in parameterizing the 
equation selected to represent the viral diffusion curve and parameterizing the equation 
selected to represent the viral diffusion curve to yield the preliminary viral diffusion curve. 

Also, in the exemplary embodiment, each dot spectrogram is mapped to the viral 
diffusion curve using fractal filtering by generating a partitioned iterated fractal system IFS 
model representative of the dot spectrogram, determining affme parameters for IFS model, 
and then mapping the dot spectrogram onto the viral diffusion curve using the IFS. Before 
the dot spectrograms is mapped to the viral diffusion curve, the dot spectrograms are 
interferometrically enhanced. After the mapping, any uncertainty in the mapped coordinates 
is compensated for using non-linear information filtering. 

In accordance with a second aspect of the invention, a method is provided for 
determining the viral load within a biological sample by analyzing an output pattern of a 
biochip to which the sample is applied. In accordance with the method, a viral diffusion 
curve associated with a therapy of interest is generated and then calibrated using at least two 
viral load measurements. Then the output pattern for the sample is mapped to coordinates on 
the calibrated viral diffusion curve using fractal filtering. The viral load is determined from 
the calibrated viral diffusion curve by interpreting the coordinates of the viral diffusion 
curve. 

Apparatus embodiments are also provided. By exploiting aspects of the invention, 
disease management decisions related to disease progression, therapy and dosage 
effectiveness may be made by tracking the coordinates on the viral diffusion curve as 
successive DNA-/RNA-based microarray samples are collected and analyzed. 
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Brief Description of the Drawings 

The features, objects, and advantages of the present invention will become more 
apparent from the detailed description set forth below when taken in conjunction with the 
drawings in which like reference characters identify correspondingly throughout and 
5 wherein: 

FIG. 1 is a flow chart illustrating an exemplary method for determining the 
effectiveness of a viral therapy in accordance with the invention. 

FIG. 2 is a flow chart illustrating an exemplary method for generating Viral 
Diffusion Curves for use with the method of Fig. 1. 
10 FIG. 3 is a flow chart illustrating an exemplary method for mapping dot 

spectrograms onto Viral Diffusion Curves using fractal filtering for use with the method of 
Fig. 1. 

FIG. 4 is a block diagram illustrating the effect of the fractal filtering of Fig. 3. 
FIG. 5 is a flow chart illustrating an exemplary method for compensating for 
1 5 uncertainty using non-linear information filtering for use with the method of Fig. 1 . 



Detailed Description of Exemplary Embodiments 



With reference to the figures, exemplary method embodiments of invention will now 
20 be described. The exemplary method will be described primarily with respect to the 
determination of changes in viral loads based upon the output patterns of a hybridized 
biochip microarray using DNA samples, but principles of the invention may be applied to 
other protein-based samples or to other types of output patterns as well. 

With reference to FIG. 1, steps will be described for generating viral diffusion curves 
25 for use is processing DNA biomicroarray output patterns to determine the effectiveness of 
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therapies imposed upon a patient providing samples for which the outputs are generated. 
Then, steps will be described for processing the specific output patterns using the VDC's. 

An underlying clinical hypothesis of the exemplary method is that antiviral treatment 
should inhibit viral replication and lower an individual' s viral load from baseline or suppress 
rising values. A stationary or rising viral load after the introduction of antiviral therapy 
indicates a lack of response to the drug(s) or the development of drug resistance. The VDC 
exploits the underlying hypothesis in part by correlating the rate of disease progression to a 
sample point value such that a change in sample point indicates progression. 

The Method 

At step 100, parameters representative of viral load studies for the therapy of interest 
are input. A preliminary viral diffusion curve is generated, at step 102, based upon the viral 
load studies. The parameters representative of the viral load studies include baseline viral 
load (BVL) set point measurements at which detection is achieved, B VL at which therapy is 
recommended and viral load markers at which dosage therapy is recommended. At step 104, 
a degree of directional causality in the preliminary viral diffusion curve is calibrated to yield 
the final viral diffusion curve. 

Steps 100 - 104 are performed off-line for setting up the VDC's, These steps need be 
performed only once for a given therapy and for a given set of baseline viral load 
measurements. Thereafter, any number of DNA biomicroarray output patterns may be 
processed using the VDC's to determine the effectiveness of the therapy. Preferably, VDC's 
are generated for an entire set of therapies that may be of interest such that, for any new 
DNA biomicroarray output pattern, the effectiveness of any of the therapies can be quickly 
determined using the set of VDC's. In general, the aforementioned steps need be repeated 
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only to update the VDC's based upon new and different baseline viral load studies or if new 
therapies of interest need to be considered. 

In the following, steps will be summarized for processing DNAbiomicroarray output 
patterns using the VDC's to determine whether any therapies of interest represented by the 
VDC ' s are effective. To determine the effectiveness of therapy at least two samples of DNA 
to be analyzed are collected from a patient, preferably taken some time apart, and 
biomicroarray patterns are generated therefrom. In other cases though, the different samples 
are collected from different patients. 

The output patterns for the DNA biomicroarray are referred to herein as dot 
spectrograms. A dot spectrogram is generated using a DNA biomicroarray for each sample 
from an N by M DNA biomicroarray. An element of the array is an "oxel": o(i,j). An 
element of the dot spectrogram is a hixel : h(i,j). The dot spectrogram is represented by cell 
amplitudes given by <B(i,j) for i: 1 to N, and j: 1 to M. 

Dot spectrograms are generated from the samples taken at different times using a 
prefabricated DNA biomicroarray at step 106. The dot spectrograms are interferometrically 
enhanced at step 108. Each dot spectrogram is then mapped to coordinates on the viral 
diffusion curves using fractal filtering at step 110. After the mapping, any uncertainty in the 
mapped coordinates is compensated for at step 112 using non-linear information filtering. 

VDC coordinates are initialized at step 114, then updated in accordance with filtered 
dot spectrograms at step 116. A degree of convergence between the mapped coordinates on 
the viral diffusion curves is then determined at step 118 and a determination is made as to 
whether the therapy of interest has been effective. The determination is based upon whether 
the degree of convergence increases from one DNA sample to another. An increase in 
degree of convergence is representative of a lack of effectiveness of the therapy of interest. 
Hence, if the degree of convergence decreases, then execution proceeds to step 120, wherein 
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a signal is output indicating that the therapy is effective. If the degree of convergence 
increases, then execution proceeds to step 122 wherein VDC temporal scale matching is 
performed. Then a determination is made at step 124 whether an effectiveness time scale has 
been exceeded. If exceeded, then a conclusion is drawn that the effectiveness of the viral 
5 therapy cannot be established even if more samples are analyzed. If not exceeded, then 
execution returns to step 106 wherein another sample taken from the same patient at a latter 
time is analyzed by repeating steps 106 through 118. 

Viral Load Studies 

! ;~ Viral load studies for therapies of interest are parameterized at step 100 as follows. 

T I The therapy of interest is selected from a predetermined list of therapies for which viral load 

"'^ studies have been performed. Measurements from viral load studies are input for therapy of 

J % interest. As noted, the viral load measurements include one or more of Baseline Viral Load 

i;«s 15 (BVL) set point measurements at which detection is achieved; BVL at which therapy is 
S recommended; and VL markers at which dosage change recommended. 

Data for the viral load measurements are obtained, for example, from drug 
qualification studies on a minimum include dosages, viral limits as well as time cycles within 
which an anti-viral drug is deemed effective. The data is typically qualified with age/weight 
20 outliers and patient history. Attribute relevant to this claim is the y x or BVL L0W which 
corresponds to the lowest detection limit shown for a therapy to be effective using 
conventional assays or any other diagnostic means. BVL^low denotes the lowest threshold 
at which viral load is achieved using a nucleotide probe. Using interferometric enhancement 
technique, BVL NP „ L0W « BVL L0W . 



25 
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Generation of Viral Diffusion Curves 

Referring now to FIG* 2, the viral diffusion curves are generated as follows. An 
equation is selected for representing the VDC at step 200. Expectation and mean response 
5 parameters are determined at step 202 for use in parameterizing the selected equation. Then 
the equation selected to represent the VDC is parameterized at step 204 to yield a numerical 
representation of the VDC. These steps will now be described in greater detail. 

These steps populate a canonical machine representation, denoted as VDC(i, r, y, k) 
(which is a special case of Fokker-Planck equation) to calibrate responses from a viral load 
1 0 detection DNA- array based hybridization biomicroarray . 

i is the index for a diagnostic condition/therapy of interest, 

r denotes the parameter vector characterizing the VDC, 

X denotes the clinical endpoints vector that indicates detectability thresholds for a 
specific DNA-hybridization array implementation 
1 5 a- correspond to the uncertainty interval estimates. 

An example of an equation selected for representing the VDC is: 

^ = div(V^(x)/>)+^ , p(x,0) = p°(x) 

ot p 

The potential W(x) : 9T -> [0, oo) is a smooth function, P > 0 is 

selected constant, and p° (x) is a probability density on SR n . 
20 Preferably, the diffusion potential of the equation and the B VL data are 

such that: 

*F(x) < c . [ BVV ^ I BVL^] 
the constant c is generally set to 
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C = [number of amplitude discretization levels] 

[ log(PCR amplification factor) * avg(oligonucleotides/oxel) * Ragging 
efficiency! * Lbinding efficiency! ] 
5 Binding efficiency is difficult to quantify analytically for a 

biomicroarray device technology. Hence, for use in the above equation, an 
estimate of the binding efficiency is preferably employed. A binding 
efficiency of 30% (.3) is appropriate, though other values may alternatively 
used. Depending upon the specific biomicroarray used, the constant c 
1 0 typically ranges between 0.000 1 to 0.5. 

Expectation and Mean Response Parameters 

The expectation (p) and mean response parameters are then 
1 5 determined at step 202 for use in parameterizing the equation selected to 
represent the VDC. The expectation and mean response values are 
determined by: 1) performing conventional PCR amplification; 2) obtaining 
calibrated viral counts from the PCR amplification; 3) determining enhanced 
and normalized hybridization amplitude mean and variance values 
20 corresponding to the calibrated viral counts; and 4) matching the enhanced 
and normalized hybridization amplitude mean and variance values. 

Two synthetic amplification techniques (in addition to PCR and any 
designer tagging) are used to achieve VL estimation above the BVL limit set 
for the exemplary embodiment of the method, namely (a) readout pre- 
25 conditioning, and (b) nonlinear interferometric enhancement. Moreover, the 
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expectation match condition implies that: 

Expectation f Logd interferometricallv enhanced imag 

-v [expression set of 

interest] 

5 Expectation [preconditioned image amplitude] v [expression set of interest 

~ 1 

Variance matching is done similarly with respect to biomicroarray 
1 0 readout. The lower bound of mean response value can be given by: 

Variance \ Log(1_ interferometricallv enhanced image]) ] w [expression set of ^est] = 
1 Variance [preconditioned image amplitude] v [expression set of ' mteTest 

15 

Using the above expression, a conservative lower bound for 
interferometric enhancement is estimated for each nucleotide expression of 
interest. Since the array fabrication device is assumed to have (i) an identical 

20 oligonucleotide density per oxel and (ii) equal length oligonucleotides, the same 
mean response amplitude can be assumed. If these two assumptions are not met 
then bounds need to be individually calculated and averaged using the above 
formula. Another assumption is that the binding efficiency is statistically 
independent of the actual oligonucleotide sequence. If this assumption does not 

25 hold for the specific device technology then the binding efficiency should be 



-16- 



provided as well for each expressed sequence of interest. So the computational 
analysis method uses the analytically derived lower bounds, as computed using 
the above equation. This is a one-time calculation only and is done offline at 
design time. 

Parameterization of the VDC 



The equation selected to represent the VDC is then parameterized using 
the expectation and mean response values at step 204 to yield a numerical 
1 0 representation of the VDC using 
y(x) = p 0 + p 1 x + p 2 x 2 



subject to constraints 



15 



x~ysm 



4co 



a 



sin cot+sit) 



with constants a, P ? y and s « 1 . 

This utility of this parameterization is established as follows. The VDC 
canonical representation is based on a variational formulation of the Fokker- 
20 Planck. The Fokker-Planck (FP) equation, or forward Kolmogorov equation, 
describes the evolution of the probability density for a stochastic process 
associated with an Ito stochastic differential equation. The exemplary method 
exploits the VDC to model a physical time-dependent phenomena in which 
randomness plays a major role. The specific variant used herein is one for 
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which the drift term is given by the gradient of a potential. For a broad class of 
potentials (that correspond to statistical variability in therapy response), a time 
discrete, iterative variational is constructed whose solutions converge to the 
solution of the Fokker-Planck equation. The time-step is governed by the 
5 Wasserstein metric on probability measures. In this formulation the dynamics 
may be regarded as a gradient flux, or a steepest descent for the free energy 
with respect to the Wasserstein metric. This parameterization draws from theory 
of stochastic differential equations: wherein a (normalized) solution to a given 
Fokker-Planck equation represents the probability density for the position (or 
5 10 velocity) of a particle whose motion is described by a corresponding Ito 

i'x I 

^ stochastic differential equation (or Langevin equation) . The drift coefficient is a 

^ gradient. The method exploits "designer conditions" on the drift and diffusion 

y coefficients so that the stationary solution of a Fokker-Planck equation satisfies 

a a variational principle. It minimizes a certain convex free energy functional 

1 5 over an appropriate admissible class of probability densities. 
Z A physical analogy is to an optimal control problem which is related to 

the heating of a probe in a kiln. The goal is to control the heating process in 
such a way that the temperature inside the probe follows a certain desired 
temperature profile. The biomolecular analogy is to seek a certain property in 
20 the parameterized VDC-namely, an exponential jump in the VDC coordinate 
position for "small linear changes in the viral count". 

This method is in contradistinction to conventional calibration strategies 
which obtain a linear or superlinear shift in quantization parameter for an 
exponential shift in actual viral count. 
25 As noted, the form of FP equation chosen is 
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3£ = div (VF(x)^)+4 > = p° W 

where the potential ¥{x) : 91" -» [0, oo ) is a smooth function, (3 > 0 is 
selected constant, and p° (x) is a probability density on W. The solution p(t,x) 
is a probability density on 5R n for almost every fixed time t. So the distribution 
p(t,x) > 0 for almost every (t,x) e (0,oo) X SR n , and 

j^ a p(t,x) dx = l 

.for almost every t, e (0,oo). 

It is reasonably assumed that hybridization array device physics for the 
DNA biomicroarray (i.e., corresponding to the potential vj/) has an 
approximately linear response to the nucleotide concentration and the response 
is monotonic with bounded drift. So, 

where the partition function Z is given by 



Z = ft-f^dx 
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In this model the basis for device physics design is that the potential 
needs to be modulated such that it grows rapidly enough for Z to be finite. This 
is not achieved by conventional methods. However, a technique which does 
achieve this result is described in co-pending U.S. Patent Application 
5 xx/xxx,xxx, filed contemporaneously herewith, entitled " Method and 
Apparatus for Analyzing Hybridized DNA Microarray Patterns Using Resonant 
Interactions Employing Quantum Expressor Functions", which is incorporated 
by reference herein. 

The probability measure p s {x)dx is the unique invariant measure for 

10 the Markov random field (MRF) fit to the empirical viral load data. 

The method exploits a special dynamical effect to design p . The 
method restricts the FP equation form above to a more restricted case: random 
walk emulating between critical equilibrium points. 

To aid in understanding this aspect of the invention, consider the 
15 diffusion form 

dp{x,t) ^\ d2 d 2 p(x,t) 
dt 2 d 2 t 

where D 2 = n a 2 

20 and a is constant. 

A specific VDC shape is parameterized by: 



y(x) = p 0 + p 1 x + p 2 x 2 
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subject to constraints 



x - y sin 



a 



sin cut +s(t) 



with constants a, P, y and 8 « 1 . 

These constants are set based upon the dynamic range expected for the 
viral load. Thus, if the viral load is expected to vary only within a factor of 10, 
the constants are set accordingly. If the viral load is expected to vary within a 
greater range, different constants are employed. The actual values of the 
constants also depend upon the particular disease. 

Where the following conditions are met 

Expectation match: E{x) = £ xf(x)dx = // 

2 

Variance: a 2 = \{x - //) f x dx 



And //(*) = 1 

— 00 

The expectation and mean response parameters for use in these 
equations are derived, as described above, from matching the enhanced and 
normalized hybridization amplitude mean and variance that correspond to 



-21- 



calibrated viral counts (via classical PCR amplification). 

A distribution represented by the above-equations then satisfies the 
following form with a prescribed probability distribution 

k 

x = y sin 

Assuming: 

y = ^ > fi > 0, constant 

and <£*(/) = ay 
such that 

a = a^(y-l)(j> + l) - s 0 a 

and the distribution controlling equation is 

/(x) = 0.5|y(x)| 

such that y(-l)<X<y(l). 

The characteristic timescale of response for this system is given by 



a 



y{x) 



sin cot +s(f) 
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r* = — arccos 

CO 



£(1/3,1/3) a-jco 

\fr r 



the successive points must show a motion with characteristic timescale. 
The VDC is designed such that sampling time falls well within characteristic 
5 time. 

As noted the actual information used for populating above the 
parameters is available from the following: Baseline viral load (BVL) set point 
measurements at which detection is achieved; BVL at which therapy is 
recommended; and VL markers at which dosage change recommended. 
1 0 The following provides an example of preclinical data that is available to to assist in 

parameterizing the VDC. 

NOTIONAL VIRAL LOAD MANAGEMENT EXAMPLE 

1 5 This is a synthetic example to illustrate how data from clinical studies 

may be used to calibrate the VDC. 

Viral Load analysis studies, using conventional assays, in HIV 
progression have shown that neither gender, age, HC V co-infection, past history 
of symptomatic HIV-1 infection, duration of HIV- 1 infection nor risk group are 

20 associated with a higher risk of increasing baseline viral load (BVL) to the 
virologic end-point. However, patients with a high (BVL) between 400-6000 
copies had a 10-fold higher risk of increasing the level of viral load than 
patients with a BVL below 1500 copies/ml. Thus, baseline viral load set point 
measurements provide an important indicator for onset of disease. 
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Initiation of antiretro viral therapy is generally recommended when the 
CD4 + T-cell count is < 600 cells/mL and the viral load level is > 6,000 
copies/mL. When the viral load is > 28000 copies/mL, initiation of therapy is 
recommended regardless of other laboratory markers and clinical status. 

Effective antiretroviral treatment may be measured by changes in 
plasma HIV RNA levels. The ideal end point for effective antiviral therapy is to 
achieve undetectable levels of virus (< 400 copies/mL). A decrease in HIV 
RNA levels of at least 0.5 log suggests effective treatment, while a return to 
pretreatment values (±0.5 log) suggests failure of drug treatment. 

When HIV RNA levels decline initially but return to pretreatment 
levels, the loss of therapy effectiveness has been associated with the presence of 
drug-resistant HIV strains. 

The therapy-specific preclinical viral load markers (such as low and 
high limits in the above example) are used to establish actual BVL boundaries 
for the VDC associated with a particular therapy. In this regard, the 
deterimination of the BVL parameters is disease specific. For example, in HIV 
methods such as RT-PCR, bDNA or NASBA are used. Other diseases use other 
assays. Typically, once the parameters of the VDC equation have been set (i.e. 
constants a, p, y and s), only two viral load markers are needed to complete the 
parameterization of the VDC. This is in contrast to previous techniques 
whereby expensive and laborious techiques are required to determine the shape 
of a viral diffusion curve. The present invention succeeds in using only two 
viral load markers in most cases by exploiting the canonical VDC described 
above which has predefined properties and which is predetermined based upon 
the particular biochip being used. 
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Calibration of Viral Diffusion Curves 

Referring again to FIG* 1, the directional causality of the VDC is 
5 calibrated at step 1 04 in the context of an NIF discussed in greater detail below. 
At least arbitrarily selected three sample points are used execute the NIF 
calibration computation. The resulting polynomial is used to extracting 
qualitative coherence properties of the system. 

The spectral [0] and temporal coherence [S] is incrementally estimated 
1 0 and computed for each mutation/oligonucleotide of interest by a NIF forward 
estimation computation (described further below). The two estimates are 
normalized and convolved to yield cross-correlation function over time. The 
shape index (i.e., curvature) of the minima is used as a measure for directional 
causality. Absence of curvature divergence is used to detect high directional 
1 5 causality in the system. 



Sample Collection 

20 Samples are collected at step 106 subject to a sample point collection 

separation amount. The separation amount for two samples is preferably within 
half a "drug effectiveness mean time" covering 2a population level wherein a 
denotes the standard deviation in period before which a effectiveness for a 
particular drug is indicated. The following are some general guidelines for 

25 sample preparation for use with the exemplary method: 
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- It is important that assay specimen requirements be strictly followed to 
avoid degradation of viral KNA. 

- A baseline should be established for each patient with two specimens 
drawn two to four weeks apart. 

5 - Patients should be monitored periodically, every three to four months 

or more frequently if therapy is changed. A viral load level that remains at 
baseline or a rising level indicates a need for change in therapy. 

- Too much significance should not be given to any one viral load result. 
Only sustained increases or decreases of 0.001- 0.01 log [conventional methods 

10 typically require a 0.3-0.5 log change] or more should be considered significant. 
Biological and technical variation of up to 0.01 log [typical conventional limit: 
0.3 log] is possible. Also, recent immunization, opportunistic infections and 
other conditions may cause transient increases in viral load levels. 

- A new baseline for each patient should be established when changing 
1 5 laboratories or methods. 



Recommendations for frequency of testing are as follows: 

20 - establish baseline: 2 measurements, 2 to 4 weeks apart 

- every 3 to 4 months or in conjunction with CD4 + T-cell counts 

- 3 to 4 weeks after initiating or changing antiretroviral therapy 

- shorter intervals as critical decisions are made. 

- measurements 2-3 weeks apart to determine a baseline measurement. 
25 - repeat every 3-6 months thereafter in conjunction with CD4 counts to 
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monitor viral load and T-cell count. 

- avoid viral load measurements for 3-4 weeks following an 
immunization or within one month of an infection. 

- a new baseline for each patient should be established when changing 
5 laboratories or methods. 

The samples are applied to a prefabricated DNA biomicroarray to 
generate one or more dot spectrograms each denoted <D(i j) for i: 1 to N, and j : 
1 to M. The first sample is referred to herein as the k=l sample, the second as 
the k=2 sample, and so on. 

10 

Interferometric Enhancement of the Dot Spectrogram 

Each dot spectrogram provided by the DNA biomicroarray is filtered at 
step 108 to yield enhanced dot spectrograms <£(k) either by performing a 

1 5 conventional Nucleic Assay Amplification or by applying preconditioning and 
normalization steps as described in the co-pending Patent Application having 
Serial Number xx/xxx,xxx, (Attorney Docket Number P67 41736), entitled 
"Method And Apparatus For Analyzing Hybridized Biochip Patterns Using 
Resonance Interactions Employing Quantum Expressor Functions". The 

20 application is incorporated by reference herein, particularly insofar as the 
descriptions of the use of preconditioning and normalization curves are 
concerned. 

Fractal Filtering 

25 
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Each enhanced dot spectrogram is then mapped to the VDC using fractal 
filtering at step 110 as shown in FIG. 3 by generating a partitioned iterated 
fractal system 302, determining affine parameters for the IFS 304 and then 
mapping the enhanced dot spectrogram onto the VDC using the IFS, step 306. 



The VDC representation models a stochastic process given by 



l 



... "N 



+- 



+ T 



.J \ 



\y-yo. 



if( x >y) e M l 1 (1) » for some 1 - i - m ' 



W{f){x,y) = - 



0, 



otherwise; 



for any (x,y)e$R 2 and f e p(9t 2 ) 
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An exemplary partitioned iterated fractal system (IFS) model for the 
system is 

where the affine parameters for the IFS transformation are given by 



T = 



V 5 io S \\J 



where the D - origin is given by (x l D ,y l D ) 9 
the R - origin is given by (x l R , y R ) 
spatial transformation matrix is given by or. 
the intensity tilting vector is given by t x 
the contrast adjuctment is given by y u 
the brightness adjustment is given by /? 



and wherein <D represents the enhanced dot spectrogram and wherein ju 
1 0 represents the calculated expectation match values 

This IFS model maps the dot spectrogram to a point on the VDC 
wherein each VDC coordinate is denoted by VDC(t,0) such that 



15 



W[0>,k] -> VDC(k ? @) 



Wherein k represents a sample. 



In the above equation, 0 represent the parameters of the IFS map. 
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Thus the output of step 110, is a set of VDC coordinates, identified as 
VDC(k, 0), with one set of coordinates for each enhanced dot spectrogram k = 
l,2...n. 

The effect of the steps of FIG. 3 is illustrated in FIG. 4 which shows a 
5 set of dot spectrograms 450, 451 and 452 and a VDC 454. As illustrated, each 
dot spectrogram is mapped to a point on the VDC. Convergence toward a 
single point on the VDC implies ineffectiveness of the viral therapy. A 
convergence test is described below. 

10 Uncertainty Compensation 

With reference to FIG. 5, any uncertainty in the coordinates VDC(k, G>) 
is compensated using Non-linear Information Filtering as follows. 
Biomicroarray dispersion coefficients, hybridization process variability values 
15 and empirical variance are determined at step 402. The biomicroarray 
dispersion coefficients, hybridization process variability values and empirical 
variance are then converted at step 404 to parameters for use in NIF. The NIF 
is then applied at step 406 to the VDC coordinates generated at step 106 of 
FIG. 1. 

20 Nonlinear information filter (NIF), is a nonlinear variant of the Extended 

Kalman Filter. A nonlinear system is considered. Linearizing the state and 
observation equations, a linear estimator which keeps track of total state 
estimates is provided. The linearized parameters and filter equations are 
expressed in information space. This gives a filter that predicts and estimates 



25 
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information about nonlinear state parameters given nonlinear observations and 
nonlinear system dynamics. 

The information Filter (IF) is essentially a Kalman Filter expressed in 
terms of measures of the amount of about the parameter of interest instead of 
tracking the states themselves, i.e., track the inverse covariance form of the 
Kalman filter. Information here is in the Fisher sense, i.e. a measure of the 
amount of information about a parameter present in the observations. 

Uncertainty bars are estimated using NIF algorithm. The parameters 
depend on biomicroarray dispersion coefficients, hybridization process 
variability and empirical variance indicated in the trial studies. 

One particular advantage of the method of the invention is that it can 
also be used to capture the dispersion from individual to individual, therapy to 
therapy etc. It is extremely useful and enabling to the method in that it can be 
apriori analytically set to a prechosen value and can be used to control the 
quality of biomicroarray output mapping to VDC coordinates. 

The biomicroarray dispersion coefficients, hybridization process 
variability values and empirical variance are determined as follows. Palm 
generator functions are used to capture stochastic variability in hybridization 
binding efficacy. This method draws upon results in stochastic integral 
geometry and geometric probability theory. 

Geometric measures are constructed to estimate and bound the 
amplitude wanderings to facilitate detection. In particular we seek a measure for 
each mutation-recognizer centered (MRC-) hixel that is invariant to local 



-31- 

degradation. Measure which can be expressed by multiple integrals of the form 
m{Z)=[f{z)dz 

5 where Z denotes the set of mutations of interest. In other words, we 

determine the function f(z) under the condition that m(z) should be invariant 
with respect to all dispersions Also, up to a constant factor, this measure is 
the only one which is invariant under a group of motions in a plane. In 
principle, we derive deterministic analytical transformations on each MRC- 
10 hixel., that map error-elliptic dispersion bound defined on SR 2 (the two 
dimension Euclidean space - i.e., oxel layout) onto measures defined on $R. The 
dispersion bound is given by 

Lo g4 (6 ( ^ r>. 

15 

Such a representation of uniqueness facilitates the rapid decimation of 
the search space. It is implemented using a filter constructed using measure- 
theoretic arguments. The transformation under consideration has its theoretical 
basis in the Palm Distribution Theory for point processes in Euclidean spaces, 
20 as well as in a new treatment in the problem of probabilistic description of 
MRC-hixel dispersion generated by a geometrical processes. Latter is reduced 
to a calculation of intensities of point processes. Recall that a point process in 
some product space E X F is a collection of random realizations of that space 
represented as {(e 1? fj), | ^ e E, fj e F }. 
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The Palm distribution, n of a translation (T J invariant, finite intensity 
point process in 9t n is defined to the conditional distribution of the process. Its 
importance is rooted in the fact that it provides a complete probabilistic 
description of a geometrical process. 

In the general form, the Palm distribution can be expressed in terms of 
a Lebesgue factorization of the form 

E P N*= ALh x n 

Where II and A completely and uniquely determine the source 
distribution P of the translation invariant point process. Also, E P N* denotes 
the first moment measure of the point process and 1^ is a probability measure. 

Thus a determination of II and A is needed which can uniquely encode 
the dispersion and amplitude wandering associated with the MRC-hixel. This is 
achieved by solving a set of equations involving Palm Distribution for each 
hybridization (i.e., mutation of interest). Each hybridization is treated as a 
manifestation of a stochastic point process in SR 2 

In order to determine II and A we have implemented the following 
measure-theoretic filter: 

Determination of A 

using integral formulae constructed using the marginal density functions 
for the point spread associated with MRC-hixel(i,j) 

The oligonucleotide density per oxel p m(hj) , PCR amplification protocol 

(ctJ, fluorescence binding efficiency (r)^ and imaging performance (w m ) 



-33- 



provide the continuous probability density function for amplitude wandering in 
the m-th MRC-hixel of interest. Let this distribution be given by 

-The method requires a preset binding dispersion limit to be provided to 



compute A. The second moment to the function 

at SNR = 0 condition is used to provide the bound. 

Determination of II 

Obtained by solving the inverse problem 
I7=@*P 

where 



where x Y and t 2 represent the normalized hybridization dispersion limits. 

The number are empirically plugged in. The values of 0.1 and 0.7 are 
appropriate for, respectively, signifying loss of 1 0% - 70% hybridization. Also 
, ® denotes the distribution of known point process. The 
forml/(l+exp( p (...))) is employed herein to represent 0 . 
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The biomicroarray dispersion coefficients, hybridization process 
variability values and empirical variance are then converted to parameters at 
step 304 for use in NIF as follows. 

The NIF is represented by: 

Predicted State = f(current state, observation model, information uncertainty, 
information model) 

Detailed equations are given below. 

In the biomicroarray context, NIF is an information-theoretic filter that 
predicts and estimates information about nonlinear state parameters (quality of 
observable) given nonlinear observations (e.g., post hybridization imaging) and 
nonlinear system dynamics (spatio-temporal hybridization degradation). The 
NIF is expressed in terms of measures of the amount of information about the 
observable (i.e., parameter of interest) instead of tracking the states themselves. 
It has been defined as the inverse covariance of the Kalman filter, where the 
information is in the Fisher sense, i.e, a measure of the amount of information 
about o r present in the observations Z k where the Fisher information matrix is 
the covariance of the score function. 

In a classical sense the biomicroarray output samples can be described 
by the nonlinear discrete-time state transition equation of the form: 



VDC(k) = f(VDC(k-l\0{k-\\k) + v(k) 
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where VDC(k-l) is the state at time instant (k-1), 
O(k-l) is the input vector (embodied by dosage and/or therapy) 
v(k) is some additive noise; corresponds to the biomicroarray dispersion 
as computed by the Palm Generator functions above. 
5 VDC(k) is the state at time k, 

f(k,.,) is the nonlinear state transition function mapping previous state 
and current input to the current state. In this case it is the fractal mapping that 
provides the VDC coordinate at time k. 
g;n The observations of the state of the system are made according to a non- 

Ijf 10 linear observation equation of the form 
W z(k) = h(VDC(k)) + w(k) 

ly where z(k) is the observation made at time k 

* VDC(k) is the state at time k, 

itf w(k) is some additive observation noise 

& 1 5 and h(. ,k) is the current non-linear observation model mapping current 

^ state to observations, i.e., sequence-by-hybridization made at time k, 

v(k) and w(k) are temporally uncorrelated and zero-mean. This is true 
for the biomicroarray in how protocol uncertainties, binding dynamics and 
hybridization degradation are unrelated and additive. The process and 
20 observation noises are uncorrelated. 

E[v(i) w T (j)] = 0 ? Vij. 



The dispersion coefficients together define the nonlinear observation 

25 model. 
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The nonlinear information prediction equation is given by 
y(k\k-T) = Y(k\k-\)f{k,VDC{k-l\k-\),u{k-l)) 
5 Y{k\k-\)^f x (k)Y-\k-l\k-\) W/(k) + Q(k)\- 1 

The nonlinear estimation equations are given by 
y{k\ k) = y(k\k-l) + i(k) 
Y{k\k) = Y(k\k-l) + l(k) 
where 

I(k) = VftJ(Ar) R' 1 (k) Vh x (k) 

i(k) = Vh T x (k) R- 1 (k) [ v(k) + Vh x (k) VDC{k jfc - 1)] 
where 

20 v(k) = z(k)-h(x(k\k-l)). 




In this method NIF helps to bind the variability in the VDC coordinate 
mapping across sample to sample so that dosage and therapy effectiveness can 
be accurately tracked. 
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The NIF is then applied to enhanced, fractal-filtered dot spectrogram at 
step 306 as follows. States being tracked correspond to post-hybridization dot 
spectrogram in this method. NIF computation as described above specifies the 
order interval estimate associated with a VDC point. It will explain and bound 
5 the variability in Viral load estimations for the same patient from laboratory to 
laboratory. 

The NIF also specifies how accurate each VDC coordinate is given the 
observation model and nucleotide set being analyzed. 

10 Convergence Testing 

Referring again to FIG. 1, once any uncertainty is compensated, the 
VDC coordinates are renormalized at step 114. The renormalized VDC 
coordinates are patient specific and therapy specific. Alternately the coordinates 
15 could be virus / nucleotide marker specific. The NIF-compensated VDC 
coordinates are renormalized to the first diagnostic sample point obtained using 
the biomicroarray. Thus a patient can be referenced to any point on the VDC. 

This renormalization step ensures that VDC properties are maintained, 
20 notwithstanding information uncertainties as indicated by the NIF correction 
terms. The approach is drawn from "renormalization-group' 1 approach used for 
dealing with problems with many scales. In general the purpose of 
renormalization is to eliminate an energy scale, length scale or any other term 
that could produce an effective interaction with arbitrary coupling constants. 
25 The strategy is to tackle the problem in steps, one step for every length scale. In 
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this method the renormalization methodology is abstracted and applied during a 
posteriori regularization to incorporate information uncertainty and sample-to- 
sample variations. 

This is in contradistinction to current viral load measurement calibration 
5 methods that either generate samples with same protocol and same assumptions 
of uncertainty or use some constant correction term. Both existing approaches 
skew the viral load readout so that measurements are actually accurate only in a 
limited "information" and "observability" context. This explains the large 
variations in readings from different laboratories and technicians for the same 

10 patient sample. 

Specifically, we include the dynamic NIF correction function to the 
gradient of the VDC at the sample point normalized in a manner such that when 
the information uncertainty is null, the correction term vanishes. As discussed 
in the above steps, the NIF correction terms is actually derived from the noise 

1 5 statistics of the microarray sample. 

<VDC(K@)> = VDC(h@) + [VNIF(YJ) k ] 

where VNIF(Y,I) k denotes the gradient of nonlinear information 
20 prediction function. Under perfect observation model this term vanishes. 

Once initialized, the VDC coordinates are then updated at step 116 
applying the IFS filter W[ ] on k+lth sample, by 



25 



VDC(k+l, ®)<r W[Biomicroarray Output, K+l]; 
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A direction convergence test is next performed at step 118 to determine 
whether the selected therapy has been effective. If convergence establishes that 
the viral load for the patient is moving in a direction representative of a lower 
viral load, then the therapy is deemed effective. The system is deemed to be 
converging toward a lower viral load if and only if: 



VDCitJ-VDCit^) 



VDC(t k ^)-VDC(t^) 



> 1 



VDC peak -VDC(t k ) 



VDC peak ~VDC{t k _ x ) 



< 1 fork>2 and j>0 



The above relationships needs to be monotonically persistent for at least 
1 0 two combinations of k and j . 

Also, date[k] - date[j] < k * characteristic time, f (in days) 



O 15 



20 



Where k captures the population variability. Typically, k < 1.2 . 

The peak VDC value is determined based on the VDC. The peak 
amplitude is an artifact of the specific parameterization to the Fokker-Planck 
equation used in deriving the VDC. It is almost always derived independent of 
the specific sample. 

In connection with step 118, a VDC Shift factor A may be specified at 
which a dosage effectiveness decision and/or disease progression decision can 
be made. The VDC shift factor is applied to estimate the VDC curvature 
traversed between two measurements. 
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If the system is deemed to be converging toward a lower viral load, an 
output signal is generated indicating that the therapy of interest is effective at 
step 120. If not, then the execution proceeds to step 122 wherein VDC scale 
matching is performed. A key assumption underlying this method is that 
movement along VDC is significant if and only if the sample points are with a 
constant multiple of temporal scale characterizing the VDC. This does not in 
any way preclude the pharmacological relevance associated with the datapoints. 
But complete pharmacological interpretation of the sample points is outside the 
scope of this method. The process is assumed to be cyclostationary or at a 
large time scale and two or more sample points have been mapped to VDC 
coordinates. The coordinates are then plugged into an analytic to estimate the 
empirical cycle time (?). This is implemented as described in the following 
sections. 

Again the empirical cycle time (?) is used to establish decision 
convergence. 

Scale Matching 

Select a forcing function of the form: 
*f / = k. p m cos<yt 

where k is a constant and m is a small odd integer (m<7). 
The phase space for this dynamical system is represented by: 
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k 



x = y sm 




a 



erf m 




sin^t 



where 



-1 



if x < -N 



erf m(x) — < 



erf(x) 
1 



if | x | < N 
x>N 



and erf m(.) denotes the error function. 

K is set to 1 ; where o< a, y, o < 1 are refer to constants. 

Let denote the cycle-time-scale for this empirical system. 

If log e (x emp /T) > 1 (in step 10) then we claim that time-scales do not 

match. 

Time Scale Testing 

Next a determination has been made as to whether an effectiveness 
timescale has been exceeded at step 124 by: 

Checking if a time step between successive sampling has exceeded T by 

determining if Time k+1 - Time k > T 

such that VDC(k+l , 0)- VDC(k, 0) < q where 

is set to 0.0001 and wherein 
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T is given by 



T* = — arccos 



' _ £(1/3,1/3) a47o 



r 



5 

B(l/3,l/3) represents the Beta function around the coordinates (1/3,1/3), 
We can actually use all B(l/2i+l,l/2i+l) for i>l and i<7. 

If Time k+1 - Time k > T then output signal at step 126 indicating that 

10 either 

no change in viral load concluded, OR 
therapy deemed ineffective, OR 
dosage deemed suboptimal. 

15 

If Time k+1 - Time k < T then process another sample by repeating all 
steps beginning with Step 4 wherein a dot spectrogram is generated for a new 
sample. 

If the effectiveness time scale has been exceeded then a signal is output 
20 indicating that no determination can be as to whether the therapy of interest is 
effective. If the time scale is not exceeded, then execution returns to step 106 
for processing another sample. If available, and the processing steps are 
repeated. 
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Alternative Implementations 

Details regarding a related implementation may be found in co-pending 
U.S. Patent Application xx/xxx,xxx, (Attorney Docket No. P67 41854), filed 
contemporaneously herewith, entitled "Method and Apparatus for Interpreting 
Hybridized Bioelectronic DNA Microarray Patterns Using Self Scaling 
Convergent Reverberant Dynamics", and incorporated by reference herein. 

The exemplary embodiments have been primarily described with 
reference to flow charts illustrating pertinent features of the embodiments. Each 
method step also represents a hardware or software component for performing 
the corresponding step. These components are also referred to herein as a 
"means for" performing the step. It should be appreciated that not all 
components of a complete implementation of a practical system are necessarily 
illustrated or described in detail. Rather, only those components necessary for a 
thorough understanding of the invention have been illustrated and described in 
detail. Actual implementations may contain more components or, depending 
upon the implementation, may contain fewer components. 

The description of the exemplary embodiments is provided to enable 
any person skilled in the art to make or use the present invention. Various 
modifications to these embodiments will be readily apparent to those skilled in 
the art and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the invention is 
not intended to be limited to the embodiments shown herein but is to be 
accorded the widest scope consistent with the principles and novel features 
disclosed herein. 
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Claims: 

WHAT IS CLAIMED IS: 

5 LA method for determining the effectiveness of a therapy by 

analyzing biochip output patterns generated from biological samples taken at 
different sampling times from a patient undergoing the therapy, said method 
||f comprising the steps of: 

W generating a viral diffusion curve associated with a therapy of interest; 

1 6 enhancing the output pattern interfermotetrically by performing a nucleic 
JU acid amplification; 

1% mapping each of the enhanced output patterns representative of 

i;5 hybridization activity to respective coordinates on the viral diffusion curve using 
fractal filtering; 

1 5 determining a degree of convergence between the mapped coordinates on 

the viral diffusion curves; and 

determining whether the therapy of interest has been effective based upon 
the degree of convergence from one sample to another. 

20 2. A method for determining the effectiveness of a therapy by analyzing 
biochip output patterns generated from biological samples taken at different 
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sampling times from a patient undergoing the therapy, said method comprising the 
steps of: 

generating a viral diffusion curve associated with a therapy of interest; 
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mapping each of the output patterns representative of hybridization activity 
to respective coordinates on the viral diffusion curve using fractal filtering; 

determining a degree of convergence between the mapped coordinates on 
the viral diffusion curves; 

determining whether the therapy of interest has been effective based upon 
the degree of convergence from one sample to another; and 

wherein the biological sample is selected from a group consisting of a 
DNA, RNA, protein, peptide-nucleic acid (PNA) and targeted nucleic amplification 
(TNA) samples. 
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Abstract of the Disclosure 

A technique is described for determining the effectiveness of medical therapy and 
dosage formulations by analyzing dot spectrograms representative of quantized 
hybridization activity in biological samples, such as DNA, RNA or other protein 
biomolecular array samples, taken at different sampling times from a patient 
undergoing the treatment. The technique directly lends itself to disease progression 
analysis based on markers such as viral load. In accordance with the technique, a 
viral diffusion curve associated with a therapy of interest is generated and each dot 
spectrogram is then mapped to the viral diffusion curve using fractal filtering to 
yield a filtered viral diffusion curve for each sample. A degree of convergence 
between the filtered viral diffusion curves is determined. Then, a determination is 
made as to whether the therapy of interest has been effective by determining 
whether the degree of convergence increases from one sample to another, with an 
increase in the degree of convergence being representative of a lack of 
effectiveness of the therapy of interest. In a specific example described herein, the 
viral diffusion curve is generated by inputting parameters representative of viral 
load studies for the therapy of interest, generating a preliminary viral diffusion 
curve based upon the viral load studies and then calibrating a degree of directional 
causality in the preliminary viral diffusion curve to yield the viral diffusion curve. 
Each dot spectrogram is mapped to the viral diffusion curve using fractal filtering 
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by generating a partitioned iterated fractal system (IFS) model representative of the 
dot spectrogram, determining affme parameters for IFS model, and then mapping 
the dot spectrogram onto the viral diffusion curve using the IFS. Before the dot 
spectrogram is mapped to the viral diffusion curve, the dot spectrogram is 
interferometrically enhanced. After the mapping, any uncertainty in the filtered 
viral diffusion curve is compensated for using non-linear information filtering. A 
method is also described for determining the viral load within a biological sample 
by analyzing a dot spectrogram generated for the sample in connection with viral 
diffusion curves associated with a therapy of interest. Thus a technique is provided 
for detecting and tracking infections such as viral infections and establishing 
clinical endpoints, based on accurate biomolecular measurements of viral DNA or 
RNA in peripheral blood. The technique also provides a computational protocol for 
leveraging a clinical marker to establish and track therapy effectiveness based on 
quantification of amplified nucleic acid i.e., DNA and RNA assays. The technique 
can potentially amplify dynamic ranges over 100X or more over conventional 
assays. The technique enables point-of-care viral load detection biosensors to 
reliably predict the likelihood of disease progression and thereby allows the patient 
to make earlier and more effective decisions about treatment. 
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Declaration and Power of Attorney For Patent Application 

English Language Declaration 



As a below named inventor, I hereby declare that: 



My residence, post office address and citizenship are as stated below next to my name, 

I believe I am the original, first and sole inventor (if only one name is listed below) or an original, 
first and joint inventor (if plural names are listed below) of the subject matter which is claimed and for 
which a patent is sought on the invention entitled 

METHOD AND APPARATUS FOR EXPONENTIALLY CONVERGENT THERAPY EFFECTIVENESS 
MONITORING USING 1)NA MTCKOARRAY BASED VIRAL LOAD MEASUREMENTS 

O the specification of which 

(check one) 

%M is is attached hereto. 

ill □ was filed on 

* Application Number 

3*3 and was amended on 

y (if applicable) 

i;g I hereby state that I have reviewed and understand the contents of the above identified specification, 

© including the claims, as amended by any amendment referred to above, 

I acknowledge the duty to disclose to the United States Patent and Trademark Office all information 
known to me to be material to patentability as defined in Title 37 , Code of Federal Regulations, 
Section 1.56. 



as United States Application No, or PCT International 



I hereby claim foreign priority benefits under Title 35, United States Code, Section 119(a)-(d) or 
Section 365(b) of any foreign application(s) for patent or inventor's certificate, or Section 365(a) of 
any PCT International application which designated at least one country other than the United States, 
listed below and have also identified below, by checking the box, any foreign application for patent or 
inventor's certificate or PCT International application having a filing date before that of the application 
on which priority is claimed. 



Prior Foreign Application(s) Priority Not Claimed 
□ 



(Number) (Country) (Day/Month/Year Filed) 

□ 

(Number) (Country) (Day/Month/Year Filed) 

■ — □ 



(Number) (Country) (Day/MontlVYear Filed) 



SDN 
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I hereby claim the benefit under 35 U.S.C. Section 119(e) of any United States provisional 
application(s) listed below; 



(Application Serial No.) 



(Filing Date) 



(Application Serial No.) 



(Filing Date) 



(Application Serial No.) 



(Filing Date) 



I hereby claim the benefit under 35 U. S. C. Section 120 of any United States application(s), or 
Section 365(c) of any PCT International application designating the United States, listed below and, 
insofar as the subject matter of each of the claims of this application is not disclosed in the prior 
United States or PCT International application in the manner provided by the first paragraph of 35 
U.S.C. Section 112, I acknowledge the duty to disclose to the United States Patent and Trademark 
Office all information known to me to be material to patentability as defined in Title 37, C. F..R., 
Section 1.56 which became available between the filing date of the prior application and the national 
or PCT International filing date of this application: 



{Application Serial No.) 



(Filing Date) 



(Status) 

(patented, pending, abandoned) 



(Application Serial No.) 



(Filing Date) 



'if*!? 



(Status) 

(patented, pending, abandoned) 



(Application Serial No.) 



(Filing Date) 



(Status) 

(patented, pending, abandoned) 



I hereby declare that all statements made herein of my own knowledge are true and that all 
statements made on information and belief are believed to be true; and further that these statements 
were made with the knowledge that willful false statements and the like so made are punishable by 
fine or imprisonment, or both, under Section 1001 of Title 18 of the United States Code and that such 
willful false statements may jeopardize the validity of the application or any patent issued thereon. 
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POWER OF ATTORNEY: As a named inventor, I hereby appoint the following attorney(s) and/or 
agent(s) to prosecute this application and transact all business in the Patent and Trademark Office 
connected therewith, (list name and registration number) 

Hank! CJ^all^Xvq^ PlLLRcgiiilraliflii-N.o^23,S22 



Send Correspondence to: Dankl c Ma,ler * 
O PRETTY, SCHROEDER & POPLAWSKT 

Jffi 444 South Flower Street, 19th Floor 

H! Los Angeles, California 90071-2909 

W 

y Direct Telephone Calls to: (name and telephone number) 

IH D aniel C Mallery, Esq ., PRETTY, SCHROEDER & POPLAWSK I, Tel.: (213) 622-7700 



Full name of sole or first inventor 






Poonam S. Gulati 






Sole orfirst inventor's signature 


Date 






02/22/99 




Residence 






5467 La Forest Drive, La Canada, CA 9 101 1 






Citizenship 






India 






Post Office Address 






Same as above. 













Full name of second inventor, If any 




Sandecp Gulati 




SecorplkventoiteL signature 


Date 




02/22/99 


Residencej\ p ~ i 1 




5467 La FQTjst Drive, La Canada, CA 91011 




Citizenship 




India 




Post Office Address 




Same as above. 







